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ABSTRACT 

We examine the rise and sudden demise of comet C/2010 XI (Elenin) on its approach 
to perihelion. Discovered inbound at 4.2 AU, this long-period comet was predicted to 
become very bright when near perihelion, at 0.48 AU on 2011 September 10. Obser¬ 
vations starting 2011 February (heliocentric distance ~3.5 AU) indeed show the comet 
to brighten by about 11 magnitudes, with most of the increase occurring inside 1 AU 
from the Sun. The peak brightness reached m./j = 6 on UT 2011 August 12.95A0.50, 
when at ~0.83 AU from the Sun. Thereafter, the comet faded even as the heliocentric 
distance continued to decrease. We find that most of the surge in brightness in mid- 
August resulted from dust particle forward-scattering, not from a sudden increase in 
the activity. A much smaller (~3 magnitudes) brightening began on UT 2011 August 
18±1 (heliocentric distance 0.74 AU), reached a maximum on UT 2011 August 30±1 
(at 0.56 AU), and reflects the true break-up of the nucleus. This second peak was 
matched by a change in the morphology from centrally condensed to diffuse. The esti¬ 
mated cross-section of the nucleus when at 1 AU inbound was ~1 km 2 , corresponding 
to an equal-area circle of radius 0.6 km. Observations were taken after the second peak 
using the Canada-France-Hawaii 3.6-m telescope to search for surviving fragments of 
the nucleus. None were found to a limiting red magnitude r' = 24.4, corresponding to 
radii <40 m (red geometric albedo = 0.04 assumed). The brightening, the progressive 
elongation of the debris cloud and the absence of a central condensation in data taken 
after UT 2011 August 30 are consistent with disintegration of the nucleus into a power 
law size distribution of fragments with index q = 3.3A0.2 combined with the action 
of radiation pressure. In such a distribution, the largest particles contain most of the 
mass while the smallest particles dominate the scattering cross-section and apparent 
brightness. We speculate about physical processes that might cause nucleus disruption 
in a comet when still 0.7 AU from the Sun. Tidal stresses and devolatilization of the 
nucleus by sublimation are both negligible at this distance. However, the torque caused 
by mass loss, even at the very low rates measured in comet Elenin, is potentially large 
enough to be responsible by driving the nucleus to rotational instability. 



Subject headings: asteroids, comets 
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1. INTRODUCTION 


Comet C/2010 XI (Elenin) (hereafter simply “Elenin”) was discovered on UT 2010 December 
10 as a ~19.5 red magnitude object at heliocentric distance rjj = 4.221 AU (Elenin et al. 2010). 


The orbit was soon determined to be that of a long-period comet (semimajor axis, eccentricity 
and inclination were -7532 AU, 1.000064 and 1.84°, respectively) with perihelion expected at q 
= 0.482 AU on UT 2011 September 10, to be followed soon after by a close approach to Earth 
(minimum geocentric distance 0.234 AU on UT 2011 October 17). Simple (but unphysical) power- 
law extrapolations of the apparent brightness from the discovery epoch gave rise to predictions that 
Elenin would become a bright naked-eye comet near and soon after perihelion. These predictions 
at first seemed to be substantiated, as ground-based observers reported rapid brightening in early 


August to peak visual magnitudes near 8 (Sostero et al. 2011). However, Elenin was seen to be 


fainter on or about UT 2011 August 17 while the morphology was reported to change from centrally 
condensed to increasingly diffuse. The comet was last recorded in unpublished observations by 
amateur astronomers as an extremely diffuse, elongated nebulosity on UT 2011 October 23, when 
outbound at r h = 1-06 AU. In due course, comet Elenin disappeared. 

Other promising long period comets (for example, C/2012 SI (ISON) and C/2011 W3 (Love- 
joy)) have received considerable observational attention (e.g. 


Sekanina & Chodas 2012 Li et al. 


2013. Knight & Battams 2014). The disintegration of these comets is relatively easily understood 


as a consequence of intense solar heating at their small perihelion distances (0.012 AU and 0.006 
AU, respectively). Comet Elenin, however, disintegrated even before reaching its relatively distant 
perihelion at 0.48 AU, under the action of incident solar fluxes >5000 times smaller than those 
experienced by ISON and Lovejoy. Elenin therefore offers a different perspective on the mecha¬ 
nism responsible for the destruction of a long period comet on its approach to perihelion. Fading 
is recognized as important for understanding the steady-state population of long-period comets 
(specifically, it is needed to reconcile (Oort (1950) dynamical model with the observed distribution 


of cometary orbital binding energies; Wiegert & Tremaine (1999)), but the mechanism responsible 
for fading is neither well-documented nor well-understood. 


Despite widespread initial excitement and the spectacular physical development of Elenin, we 
are aware of only one relevant scientific publication in the refereed journals (Korsun et al. 2012). 
In the present paper, we present calibrated photometric measurements from a range of space-based 
and ground-based telescopes, documenting the demise of comet Elenin. We use these measurements 
to interpret the mechanism behind its disappearance. 
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OBSERVATIONS 


We obtained comet observations using a mixture of ground-based and space-based telescopes 
between 2011 February and September in the inbound-lag (see Table[l]). Ground-based photometric 
monitoring observations of the approach to perihelion were obtained between UT 2011 February 
and 2011 June using Berkeley’s Katzman Automatic Imaging Telescope (KAIT). For part of the 
period of interest, the solar elongation of comet Elenin as seen from Earth fell below 30°, effectively 
precluding ground-based observations. However, the comet entered the fields of view of cameras 
onboard the Sun-observing spacecraft STEREO A and B between 2011 May 16 and September 22. 
These spacecraft move in Earth-like orbits but are separated from the Earth in longitude by large 


angles (Howard et al. 2008), providing data from non-terrestrial perspectives and complementing 


the observations from Earth. In addition to their regular Sun-pointed observations, a set of targeted 
measurements was obtained using the STEREO B spacecraft. After the comet reached its peak- 
brightness in August and became faint in subsequent month prior to perihelion, the Canada-France- 
Hawaii Telescope was used to conduct a sensitive search for nucleus fragments in October. Lastly, 
archival measurements from the New Technology Telescope offered views of the comet at two epochs 
bracketing peak brightness at the end of July and the end of August. Figures 0 and <© show 
the comet geometry relative to the ground- and space-based telescopes between 2011 February and 
September as the comet approached the Sun. 


2.1. Katzman Automatic Imaging Telescope (KAIT) 


The KAIT is a 0.76 meter telescope located at Lick Observatory, California, at an altitude of 
1283 m (Filippenk o et al.|200 1). The telescope is mostly used for supernova survey observations. We 
used KAIT to monitor the brightness evolution of Elenin between UT 2011 February 10 and June 
3. The KAIT CCD has dimensions 512x512 pixels (24 fin i per pixel) with a field of view 6'.8x6'.8 
(at 0".8 pixel -1 ). The telescope tracks at sidereal rates without guiding, limiting the maximum 
usable exposures. Comet Elenin was observed typically in a sequence of ten images, each of 30 
s duration, for a total integration of 300 s per night in the R-band. The comet trailing does not 
substantially affect the photometry. The comet moved with the maximum speed ~54"hr~ 1 , leading 
to negligible image trailing (~ 0.5" in 30 seconds). We used a photometry aperture 7.5 pixels (6") 
in radius for the object and an annulus with inner and outer radii 9 and 12 pixels (7".2 to 9".6), 
respectively, to define the sky background. Photometric calibration was secured using Sun-like 
field stars, identified with the Aladin Sky Atlas and with R-magnitudes taken from the USNO-A2.0 
catalog. The magnitude uncertainty, typically ±0.1 magnitudes, was estimated from the scatter of 
10 individual measurements of field stars. 
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2 . 2 . STEREO Heliospheric Imagers 


The comet was detected during both regular and targeted observations by the Heliospheric 
Imagers (HI) on board the twin solar spacecraft STEREO A and B (Eyles et al. 2009). On each 
spacecraft, HI is equipped with two cameras, HI-1 and HI-2, having small (20°) and large (70°) 
angular fields of view, respectively. The camera centers point away from the Sun by 14° (HI-1) and 
53.7° (HI-2). The spectral passbands of the HI-1 cameras cover wavelengths from 630 to 730 nm, 
similar to the passband of the Johnson R-band filter. However, the pre-launch calibration indicates 
about 25% leak from the wavelengths at 300 - 450 nm (Bewsher et al. 2010). This may result a 
slight departure from the photometry from R-band. 

The STEREO HI data are publicly available and were downloaded from the UK Solar System 
Data Center. We retrieved Level 1.0 data and proceeded to remove the background static coronal 
brightness using the same median procedure that we developed to study the asteroid Phaethon 
(Li &; Jewitt; 2013). The comet should be in the field of view of the STEREO A HI-1 camera 
between May 16 and June 20, but we did not detect it because the comet brightness was below the 
detectable threshold of the camera during this time period. It was estimated that the comet had 
the brightest magnitude ~15, while the camera threshold is about 12. The comet was observed, 
however, during the special maneuver of STEREO B. The wide-held camera, HI-2, recorded the 
comet from August 01 to 05 while the narrow held camera, HI-1, recorded it from August 06 to 
12 (designated HI-lBs). However, we found that the extremely poor angular resolution of the 
HI-2 camera rendered these data susceptible to excessive contamination by background coma and 
held objects. As a result, we elected to not include HI-2 data in the present study. During the 
regular observations, the comet entered the HI-l/STEREO B held of view between August 14 and 
September 04 (designated HI-1B); and HI-l/STEREO A between September 01 to 22 (designated 
HI-1A). Figure ^ shows the comet appearance in the camera on STEREO B. 

To increase the signal-to-noise ratio and remove background coronal structures and stars near 
the comet, we extracted the comet photometry from median images shifted to be centered on the 
comet. During the targeted observations by STEREO B, the comet was observed for an hour a day. 
We obtained hourly median images of the comet and stars between August 01 and 12 (HI-lBs). 
In regular STEREO B and A observations, we made median images every 0.5 days from August 
14 to September 2 with STEREO B (HI-1B), and every 2.4 days with STEREO A (HI-1A). The 
longer median durations were used to obtain better rejection of background stars as the comet was 
projected against the Milky Way as seen from the latter spacecraft. 

Equivalent median images were calculated for standard stars in the same way as for the comet. 
Field stars for photometric calibration were selected to be spectrally similar to the Sun and spatially 
close to the comet in the images. For the HI-1 camera, stars can be found in the SECCHI star 
catalog through an IDL routine scc_get_stars which is the part of SolarSoftWare (SSW) IDL package 
(Freeland & Handy 1998). The catalog contains more than 35,000 stars brighter than V = 10. It 
provides the V-magnitude of the stars; we converted from V to R by subtracting 0.36 (the color of 
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the Sun) from the measurements. 

The STEREO HI-1 camera has a pixel size of 70". We examined the cometary radial profile, 
and decided to use a circular aperture 3 pixels (210") in radius, with the sky background determined 
in a concentric annulus having inner and outer radii of 6 and 9 pixels (7.0 7 and 10.5^ for photom¬ 
etry. The sky annulus may be weakly contaminated by the comet tail but this contamination is 
preferable to the greatly enlarged uncertainties of measurement which are incurred by the use of 
larger annuli. The photometric uncertainties were estimated based on the field star measurements. 
These uncertainties depend on the heavily under-sampled point-spread functions of the data and 
on the residual background brightness. 


2.3. The ESO New Technology Telescope 


We searched for Elenin observations in the archive data organized by the Canadian Astronomy 
Data Center. ESO’s New Technology Telescope (NTT), a 3.6m telescope, observed the comet on 
UT 2011 July 23 and August 28 (principal investigator Olivier Hainaut). The comet was observed 


with the ESO Faint Object Spectrograph and Camera (EFOSC) (Buzzoni et al. 1984; Snodgrass 


et al. 2008) on NTT in the R-band. The EFOSC CCD field of view is 4.TX4.T, while the image 


size is 2060x2060 pixels. The images are binned 2x2 pixels in the archival data giving an effective 
pixel scale 0.24" (pixel)^ 1 . The telescope was tracking the comet. Figure [ 4 ] shows the comet images 
from each date, revealing a dramatic change in the morphology between these two observations. 


Because of the dramatic changes of the comet morphology, we used different photometry 
apertures for the two dates. The comet aperture was 140 pixels (33.6"), and the sky annulus was 
150 to 200 pixel for the July data. The aperture was 350 pixels (84.0"), and the sky annulus 
was 350 to 400 pixels for the August data. We used all bright field stars in the field of view to 
calibrate the photometry. The star photometry was extracted using apertures 15 pixels (3.6") in 
radius with a concentric sky annulus from 15 to 20 pixels. Three images of the comet were taken 
on July 23 with exposure time 30 seconds within 3 minutes (UT 23:20- 23:23), and ten images 
were taken on August 28 within 11 minutes (UT 23:20-23:31), and exposure times ranged from 
3 to 80 seconds. Assuming that the comet brightness did not change dramatically on such short 
timescales, we calculated average magnitudes of the comet observed in the July 23 and August 
28. The uncertainties are estimated from the scatter of the repeated measurements. In the July 
23 data, the coma noticably over-spilled the photometry aperture. However, experiments with 
different apertures indicate that the derived magnitude is within a few xO.l magnitudes of the 
“total” magnitude that would be obtained in an aperture of infinite radius. The measurement 
on August 28 suffers more from aperture overspill, and should be regarded as a lower-limit to 
the brightness. The two data points from NTT are plotted in Figure ([5]) along with apparent 
magnitudes from KAIT and STEREO A and B. 










2.4. Canada-France-Hawaii Telescope 


We used the 3.6 meter diameter Canada-France-Hawaii Telescope (CFHT) atop Mauna Kea, 
Hawaii, to examine Elenin on UT 2011 October 22. The MegaCam prime focus imager was used 
with the r’ filter to obtain a sequence of images at the expected position of Elenin. MegaCam 
houses 36 CCDs each containing 2048x4096 pixels of 0.187" angular size, giving a total field of 
view of Px 1°. The CFHT was tracked at non-sidereal rates during the observations so as to follow 
the motion of the comet against the stars. Simultaneous guiding and non-sidereal tracking was not 
possible and therefore we limited the exposure duration to 120 s in order to minimize image smear 
due to open-loop tracking errors. The observing conditions were excellent, with 0.6 W FWHM seeing 
and photometric skies. The observations were taken in queue-scheduled mode. 

We obtained 22 images on UT 2011 Oct 22. The CCDs within each image were combined 
using the SWarp software and compared visually on a computer. No candidate objects showing the 
expected angular motion of Elenin were found, down to a limiting magnitude estimated to be r' = 
24.4. At the expected position of Elenin, we did not detect a nucleus. 


2.5. Other Observations 


Between UT June 25.93 and August 23.39 2011, a number of amateur observers reported 


sightings of the comet (Sostero et al. 2011). Visual magnitudes of active comets, especially when 


reported by different observers using different instruments, are very difficult to interpret and so 
we have not used them in the present analysis. However, the amateur observations do provide a 
qualitative reference for the comet brightness between May and mid-July when data are available 
from neither KAIT nor STEREO. They show that the apparent brightness increased gradually 
during this time period, consistent with the slow brightening evident in earlier KAIT data, and 
reaching a peak on mid-August. 

On July 30.153 and 30.249, M. Drahus and B. Yang detected HCN emission in comet Elenin 
using the James Clerk Maxwell Telescope (Drahus et al. 2011). Assuming an isotropic production of 
gas at velocity 0.5 km s _1 and a Boltzmann distribution of energy levels at 50 K, the derived HCN 
production rate is 1.5 x 10 25 molecule s -1 . They noted that the HCN production is comparable 
to the mean level measured in comet 103P/Hartley at the same (1.07 AU) heliocentric distance in 


late 2010 (Drahus et al. 2011). At this time the estimated visual magnitude was ~9.7 (Gonzalez 


et al. (2011)). During the same period, the comet’s trajectory was nearly on the same plane as 
the Earth’s orbit and its minimum distance to STEREO B was only 0.05 AU. The spacecraft was 
literally within reach of the comet tail; the PLASTIC instrument on STEREO B directly detected 
suprathermal pickup H + and He + and singly charged water group ions produced by outgassing 


from the comet (Kucharek et al. 2013). While serving to confirm that the comet was a significant 


source of water at the end of July, it is not practical to use the plasma detections alone to estimate 
the production rate. 
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The comet was too close to the Sun to be seen from the Earth in 2011 September. A search on 
UT October 9.5 and 10.6, using the Faulks Telescope North having a field of view 10'xlO', met with 
no success (Sekanina & Guido 2011). On UT Oct. 21.38 and 21.48, the same observers reported a 
detection of a large and diffuse cloud 4'.3 from the predicted comet position in the east-northeast 
direction (position angle 77°). This observation was made with the GRAS 0.1-m f/5 APO refractor 
at the Mayhill station in New Mexico (field-of-view 3.9°x2.6°; scale 3.5" (pixel) -1 ). The large blob 
moved with the predicted comet motion in the sky. On Oct. 23.4. with the same telescope, [Sostero 


et al. (2011) confirmed the “extremely faint and diffused blob” with an extended size l°.5xl0 / .0. 


This blob appears to be the expanded and diluted remnant of the post-break-up dust cloud imaged 
at NTT on UT 2011 August 28 (Figure [4]) . 


On November 8.12, the amateur astronomer Bernard Hausler, in Germany, took images at 
the predicted location of comet Elenin. He used a 30-crn Schmidt-Cassegrain with a field of view 
29 / .7x20 / .0 and the pixel scale 1".63. The comet was not detected. 


3. RESULTS 

Our observations make it clear that Elenin did not conform to the simple power-law brightening 
model used to make optimistic predictions in 2010 and 2011. In this section, we examine the 
brightness variations and attempt to use them to diagnose physical processes occurring at the 
comet. 


3.1. The Aperture and Distance Corrections 


The apparent magnitude of the comet is shown as a function of time in Figure [5] The mag¬ 
nitudes were determined using different photometry apertures and with the comet at different 
distances from the telescopes. Therefore, the apparent magnitudes are strictly not comparable, 
since they sample different volumes of space around the nucleus and different amounts of encir¬ 
cled dust. To take account of this, we first make an “aperture correction”, in which the apparent 
magnitudes are scaled to the magnitudes which would be observed if each telescope could sample 
a fixed volume (represented by a sphere of fixed radius) around the nucleus. For this purpose, we 
arbitrarily scale the photometry to estimate the magnitude that would have been obtained if the 
comet were observed using a 100" aperture from a distance of 1 AU (corresponding to a linear 
radius of 73,000 km). 


To make the aperture correction, we need to know how the coma surface brightness, E(^>), 
varies as a function of the photometry aperture radius, qb. Isotropic coma expansion at constant 
speed would give E(0) oc (j)~ s with s = 1.0, while acceleration of the dust under a constant radiation 
pressure gives a steeper gradient, s = 1.5 (Jewitt & Meech 1987). In the former case, the encircled 
dust coma cross-section (and, therefore, the brightness) grows with radius as f 2nr<j)T,((j))d<j) oc (j), 












and in the latter as oc q i 1 / 2 . The resulting aperture correction is given by 


Am 


-2.5(2 


s) logio 



(1) 


where </> is the photometry aperture radius in arcseconds, A is the observer to comet distance in 
AU and s / 2. 


The surface brightness profile of the coma was measured by Korsun et al. (2012) on UT 2011 
March 28. They found a power-law relation with s = 1.56A0.01, close to the value expected for a 
radiation-pressure swept coma. Unfortunately, the coma was not well resolved in the KAIT data 
owing to the surface brightness sensitivity being limited by the short exposures and the brightness 
of the sky on Mount Hamilton. Surface brightness measurements in STEREO data are also limited 
by the large pixel size (70" pixel -1 ) and the bright, near-Sun background of coronal and scattered 
photospheric light. We assume that the KAIT and STEREO observations of Elenin can be well- 
represented by s = 1.5, although we strictly possess evidence validating this assumption only around 
the time of the Korsun et al. (2012) observation. The coma is fully resolved in the NTT data, where 
measurements within circular, concentric apertures show that the inner surface brightness on UT 
2011 July 23 is best described by s = 1.0. The NTT data from UT 2011 August 28 are insufficient 
to permit a meaningful determination of s, owing to the low surface brightness of the extended 
debris trail on this date. To be consistent with July data by NTT, we take s = 1.0 for the August 
NTT comet image. 


Evidence that the aperture correction is at least approximately correct is provided by Figure 
[6j in which data from different apertures and different telescopes are seen to overlap within the 
uncertainties of measurement. The KAIT, HI 1A, IB, IBs and NTT data are all broadly consistent 
in the figure even though they employed quite different photometric apertures and were taken at 
different telescope-to-comet-distances. A grave error in the aperture correction would cause a mis¬ 
match between datasets taken from different observatories. Even after the aperture correction, an 
error is incurred because distant parts of the comet are excluded by the use of a finite aperture 
including the sky contamination by the coma. This error is modest (a few tenths of a magnitude or 
less) until the middle of August, when the morphology changes from centrally-condensed to diffuse. 
A compensating factor is that the post-peak photometry used very large apertures (210" in radius) 
but it is still true that the post-peak measurements should properly be regarded as setting a lower 
limit to the integrated brightness of the whole coma. The use of still larger apertures is precluded 
by the increased sky noise introduced by the bright background of coronal and scattered light, as 
well as by imperfectly removed field objects. 


We used the aperture correction and the inverse square law to calculate the absolute magnitude 
of the comet at a given phase angle (i.e. the magnitude corrected to ru= A = 1 AU) from 


m{ 1,1, a) = m(rn, A, a) — Am — 5 log 10 (rf/A) 


(2) 
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where m(rn, A, a) is the apparent magnitude and Am is from Equation (jT]) . The aperture and 
distance-corrected magnitudes of the comet are plotted in Figure Q. Variations in the photometry 
shown there result from a combination of the effects of real changes in the scattering cross-section 
of the dust within the scaled photometry aperture and changes due to the angle-dependence of the 
efficiency of the dust scattering (the so-called “phase function”). 


3.2. The Phase Function 


As the comet approached the Sun, it entered a strongly forward-scattering geometry, reaching 
the maximum phase angle a = 172.6° as observed from STEREO B (Figure[2]). While Sun-grazing 
comets are often observed by LASCO/SOHO at moderate to large phase angles (Grynko et al. 2004), 
ground-based telescopes are rarely able to point near the Sun and so large phase angle observations 
are rare. The current record for maximum phase angle is held by P/2003 T12 observed at a = 
177.6° (Hui 2013). Observations of comets at large phase angles show that the dust is strongly 
forward-scattering ( |Ney k Merrill 1976; Ney 1982; Marcus||2007 ). 


We first fitted a Gaussian function to the portion of the aperture- and distance-corrected light 
curve near the photometric maximum, between UT August 06.3 and 18.5. The best fit gives the 
date of peak brightness as UT 2011 August 12.95. The uncertainty of the peak date is about 0.5 
days estimated by a visual examination. The maximum phase angle of the comet, a = 172°.6, 
occurred on UT 2011 August 13 as viewed from STEREO B. This coincidence between the dates of 
maximum brightness and maximum phase angle is a strong evidence that the main brightness peak 
in Figures ([5]), ([6]) and ([T]) is caused by the cometary dust phase function and the forward-scattering 


observational geometry, not by a true outburst in the comet. Sekanina k Guido (2011) inferred 


breakup on UT 2011 August 16±4, consistent with our estimate. 

The comet brightened by ~3 magnitudes due to the phase angle effect (Figure [7]). However, 
the light curve is not completely symmetric about the maximum phase angle, presumably because 
of the sudden release of debris caused by the breakup of the nucleus (as witnessed by the dramatic 
change in the morphology of the comet near this time - see Figure Q). To derive a phase function, 
we use the portion of the photometry data taken by KAIT, HI-IBs and NTT before the maximum 
phase angle of the comet. This is the time period in which we can reasonably assume that there 
was no outburst based on the observations. 


Following Marcus (2007), we adopt a compound Henyey-Greenstein (HG) function (Henyey k 
Greensteinji941) to represent the comet dust. The HG function has no physical basis, but is useful 
to provide an empirical fit to the phase functions of dust: 


m = 


A(l-g 2 ) 


47t (1 + g 2 — 2g cos @) 3 / 2 


(3) 
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where 4> is the normalized flux; 9 = 1 — a is the scattering angle (the deviation of the ray from the 
forward direction), A is the dust grain albedo, and — 1 < g < +1 is the asymmetry factor, defined 
by g = Jq 2tt&( 9) cos 9 sin 9d9. Physically, g is the scattered intensity weighted by cos 9 averaged 
over the entire solid angle. The HG function is normalized such that the integral over 47r steradians 
is unity: f^~[j^ 3>(0) sin 9d9]d(j) = 1. Large positive (negative) values of g represent strong forward 
(back) scattering. Isotropic scattering has g = 0. We fitted the aperture and distance-corrected 
flux (lo~ 0 - 4xm ( 1,1 > Q )) observed by KAIT, HI-lBs and NTT before the comet phase angle reached 
the maximum to a combined forward and back scattering HG function 


fHc{ot) = Ff®(gf) + F b &(g b ). 


( 4 ) 


In Equation Q, four free parameters are to be determined: g/ > 0 and g b < 0 are the forward 
and backward scattering ^-parameters; Ff and F b are scaling constants to the distance-corrected 
flux for the forward and back scattering components. A y; 2 fit yields values for gj = 0.926 and 
g b = —0.584; Ff = 8.075 x 10“ 4 and F b = 2.850 x 10 -4 with A=0.1. The large gj is a result of 
forward scattering by particles much larger than the wavelength of observation. For comparison, 
(mostly sub-micron) dust in the interstellar medium has a more modest gf ~ 0.5 to 0.6 (Draine 


2003). The fitted HG function is plotted in Fig (|8j) Marcus (2007) re-normalized the compound 
Henyey-Greenstein function at the phase angle a = 90° as the dust scattering function. The 
“partitioning coefficient” k in his work is equal to the ratio Ff /(Ff + F b ) in the current work. 
We obtain k = 0.74. This is smaller than the value 0.95 used by Marcus (2007) indicating that 
comets do not share a single phase function, presumably because their dust properties differ. Based 
on the best-fit HG function, the brightness of the comet that would have been observed in pure 
forward-scattering (a = 180°) is 120 times that of the backscattered brightness (a = 0°). Our best 
estimate of the pre-perihelion absolute magnitude (i.e. corrected to rjj = A = 1 AU and to a = 
0°) is 771^(1,1,0) = 11.7. 

Figure © shows how the comet brightness varied with the phase angle, this time plotted so 
as to distinguish the pre-phase-angle-maximum and post-phase-angle-maximum legs of the orbit. 
Post-phase-angle-maximum brightening relative to the phase curve is evident, with a local maximum 
near a ~ 140°. 

The relation between the absolute magnitude and the effective scattering cross section of the 


dust, C e , (in km 2 ) is by Harris (1996) 


^ = 1.5 x !0 6 10 -(o. 4mv (i,i,o)) (5) 

Pv 

where pv is the geometric albedo, which we take to be 0.04 (representative of the albedos measured 
for the nuclei of comets). 

The absolute magnitude (i.e. with both the aperture correction and the distance correction 
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taken into account) is shown as a function of time in Figure (10). The absolute magnitude remained 
unchanged at 771r(1, 1, 0) ~ 11.7±0.5 from UT 2011 February through August (40 < DOY < 230), 
in data from the KAIT, NTT and HI-lBs. Substituting mjj( 1,1,0) = 11.7 and including a small 
correction for the color (we assume V-R = 0.36), we find C e = 570±300 km 2 


In Figure (10) the main forward-scattering peak centered on August 13 is gone, but a separate 


and later peak appears at the end of August. Figure is an enlargement of Figure ( |10[ ) showing 
the STEREO and NTT data recording this peak. We again fitted this second peak with a Gaussian 
profile, finding a center on UT 2011 August 30.14 at mn(l, 1,0) = 9.8 ±0.5 and a full width at half 
maximum ~12.2 days. This second peak appears unrelated to the phase angle and, instead, we 
interpret it as a real increase in the scattering cross-section of Elenin. Substituting into Equation 
([ 5 ]), we obtain C e = 4300±900 km 2 at the maximum of this second peak. The rise towards maximum 
brightness begins near DOY 230±1 (UT 2011 August 18±1), which we take as the start of the 
disintegration of the nucleus of Elenin. 


3.3. Gas Production 


Spectral detections of CN and C3 molecules on UT 2011 March 28 (th= 2.92 AU) gave 


production rates Qcn = 1.4xl0 24 s 1 and Qc 3 = 4.2xl0 21i s x , respectively (Korsun et al. 2012). 
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A very modest production rate of hydrogen cyanide, Qhcn = 1.5xl0 25 s~ 4 , was determined from 
pre-perihelion submillimeter spectroscopy on UT 2011 July 30 at rjj= 1-07 AU (Drahus et al. 2011). 
Both CN and HCN are trace species in the comae of comets. We scaled their production rates to 


that of water using the nominal ratio QohIQcn = 320 (A’Hearn et al. 1995), finding Qh 2 o ~ 
4.5xl0 26 s” 1 (12 kg s” 1 ) at r H = 2.92 AU and Qh 2 o ~ 4.8xl0 27 s" 1 (128 kg s” 1 ) at r H = 1.07 


AU. Lovell et al. (2011) failed to detect hydroxyl (OH) emission close to perihelion (r# = 0.49 AU) 
on UT 2011 September 07, setting an upper limit to the production rate of water Qh 2 o < 3xl0 27 
s~ 4 (80 kg s~ 4 ). These modest production rates are consistent with sublimation from very limited 
areas of exposed ice, as we next discuss. On 2011 June 28 {rjj = 1.6 AU), Schleicher (private 
communication, 2015 February 6) obtained the Qoh = 10 27 ' 8 s _1 with a ±40% uncertainty, and 
Qcn = lO 25 04 s _ i. 


The resulting ratio, QohIQcn = 575±230, is consistent with the A’Hearn et 


al. (1995) value of 320 that we used to estimate the water production rate. 


We solved the energy balance equation for water ice, including energy gained and lost by 
radiation and energy consumed as latent heat in the sublimation of water. This equation is written 


= x [ e(jT 4 + L(T)F S ] (6) 

r H 

in which F & = 1360 W m 2 is the Solar constant, A is the Bond albedo, e is the emissivity, a is 
the Stefan-Boltzmann constant, L(T) is the latent heat of sublimation for water ice and F s is the 
equilibrium sublimation flux (kg m~ 2 s _1 ). The term on the left hand side describes the absorbed 
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solar power while the two terms on the right describe radiation from the surface into space and 
power used to sublimate ice at rate F s . Another, much smaller term representing conduction into 
the interior has been neglected. Dimensionless parameter y represents the ratio of the absorbing 
area on the nucleus to the area from which absorbed heat is radiated. Limiting values range from 
X = 1 (a subsolar ice patch on a non-rotating nucleus) to x = 4 (a spherical, isothermal nucleus 
in which the Sun’s heat is absorbed on nr 2 and radiated from Anr 2 ). However, measurements of 
thermal radiation from cometary nuclei show that night-side radiation is negligible, as a result of 
the low diffusivity of the surface layers (Fernandez et al. 2013). Likewise, in-situ imaging from 
spacecraft shows that the bulk of the outgassing from active comets occurs on the hot day side. 
For a spherical nucleus, day-side only emission corresponds to y = 2. Accordingly, we use % = 1 
and x = 2 to bracket the highest and lowest plausible temperatures and sublimation rates on the 
nucleus, respectively. For the other parameters we take A = 0.04, e = 0.9, a = 5.67xl0 -8 W m~ 2 
K -4 , while L(T) is obtained from Washburn (1926). 


The effective area of sublimating ice needed to supply gas at rate dM/dt is given by 


A, = 


dM / dt 


( 7 ) 


Values of A s computed from the gas emission observations and Equations ([6]) and ([T]) are listed 
in Table ([2]) and plotted against heliocentric distance in Figure (12). In the Figure, red and blue 
points distinguish sublimating areas computed using the high (y = 1) and low (y = 2) temperature 
approximations, respectively. Solutions between the red and blue lines are allowed by the data. It 
is evident that the sublimating area decreased as Elenin approached the Sun, particularly inside 
1 AU heliocentric distance. This is different from the trend expected of a long-lived, un-evolving 
source (which would plot as a horizontal line in the Figure). Beyond r# ~ 1 AU the gas production 
rates are compatible with sublimation from an area A s ~ 1 km 2 . 


3.4. The Nucleus 


The sublimating area, A s , provides a useful but imperfect estimate of the nucleus size. With 
A s = 1 km 2 , we compute the effective sublimation radius r s = ( A s /t:) 4 / 2 ~ 0.6 km. The estimate 
is imperfect because the nucleus could be larger than this but with a smaller fraction of its surface 
in sublimation. It is also possible that some fraction of the gas is produced by sublimation from 
dust in the coma, rather than from the nucleus. 


Non-detections in the CFHT data set upper limits to the surviving nucleus fragment sizes, 
through Equation ([5]). On UT 2011 October 22, the heliocentric and geocentric distances were 
rn= 1.055 AU, A = 0.246 AU, respectively, and the phase angle was a = 69.3°. For an object 
with solar colors, the Sloan r' magnitude is related to the Johnson V magnitude by r' = my — 0.16 
(Smith et al. 2002), giving a limiting magnitude from the CFHT data of my = 24.56. The phase 
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function of possible nucleus fragments is unknown, but measurements of other cometary nuclei are 
broadly compatible with 0.04±0.01 magnitudes per degree (Lamy et al. 2004). After correcting 
for the distances and phase angle, we find a limit to the absolute magnitude of surviving nucleus 
fragments of rriy( 1,1,0) > 24.72±0.7. The large uncertainty on my( 1,1, 0) reflects mainly the effect 
of the unknown phase function. We substitute into Equation ([5]) assuming py = 0.04, as is typical 
of the nuclei of comets (Lamy et al. 2004), to find an upper limit to the equivalent circular area 
radius of any surviving nucleus fragment as r max = (C e /i r) 1//2 < 40 m. Larger fragments, unless of 
much lower albedo (or perhaps of much steeper phase function), would have been detected in the 
CFHT images. 


4. DISCUSSION 
4.1. Nature of the Breakup 

We have obtained the following constraints on the breakup of the nucleus of Elenin: 

• A crude estimate of the radius of the nucleus before breakup is r n ~ 0.6 km, based on the 
measured gas production rates. 

• The largest bodies remaining after breakup had radii r max < 40 m. 

• The break up event started approximately on UT 2011 August 18±1 at th = 0.7 AU. 

• The peak cross-section of the debris produced by the breakup, as determined from the peak 
centered on August 30 at 7tjr(1, 1, 0) = 9.5T0.2 and Equation ([H]), was C e = 4300 ±900 km 2 . 

We obtain a constraint on the size distribution of the post-fragmentation particles as follows. 

We suppose that the nucleus fragments into a differential power-law size distribution in which 
n(a)da = Ta~ q da. Here, T and q are constants of the distribution, and the particles span the size 
range from a m i n to a max ■ The combined cross-section of such a power-law distribution is 


C e 


uTa 2 q da 


( 8 ) 


In a breakup, the total mass of the fragments must be equal to the initial mass of the nucleus. We 
write 


‘i 3 

p n r n 



47tT/3 

3 


a 3 q da 


( 9 ) 


where p n and r n are the nucleus density and radius, respectively, and p is the density of the ejected 
particles. We combine Equations ([8]) and Q to eliminate T, obtaining 
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3 pCe fc 

P n r n = 


hnax a 3 ~ q da 


77 


r 

Ja 


' max a 2 ~ q da 


( 10 ) 


Particles smaller than ~10 - ' m are inefficient scatterers of optical photons and will contribute 
negligibly to the measured cross-section, C e , while objects larger than 40 m should have been 
detected in the post-outburst CFHT images, but were not. Accordingly, we set a m i n = ICE' m and 
o-max = 40 m, respectively. We assume that p = p n , and then, given a max 3> a m j n , we approximate 


the solution to Equation (10) by 


r 


3 

n 


Ce f3~q\ 

77 \4-qJ 


4—q 

O’max 

3-q 


( 11 ) 


provided q ^ 3,4. We solved Equation HU by Newton-Raphson iteration using the values of C e 
and r n found earlier, to find q = 3.24. By substitution into Equation ([8]), we obtain T = 3.7xl0 7 . 
The derived value of q is relatively insensitive to the assumed input parameters. For example, 
changing r n from 0.3 km to 1.2 km changes q from 3.36 to 3.06. Changing a max , a m i n or C e even 
by an order of magnitude has a similar or smaller effect. While we make no attempt to define a 
formal statistical uncertainty, from a range of experiments we are confident that distributions with 
q = 3.3T0.2 encompass the likely range of input parameter uncertainties. 

In a q = 3.3 distribution, the mass is carried by the largest particles while the cross-section is 
dominated by small particles. To see this, we calculate ai / 2 , the particle radius below which the 
integrated cross-section is half the total cross-section, from 


1 

2 


f 

Jc 


11/2 a 2 q da 


L 


max a 2-q da ■ 
U"min 


( 12 ) 


With q = 3.3 we solve Equation (12) to find aq /2 = 1 jum. Half the cross-section, and presumably 
half the scattered intensity, is carried by particles smaller than 1 pm. However, these sub-micron 
particles carry only ~5xl0 -6 of the total mass. 


Small particles, in addition to having a strong effect on the scattered intensity, can be strongly 
affected by solar radiation pressure. Particles accelerated by a constant radiation pressure travel a 
distance from their source given by 


Sx = ^/dg Q 6t 2 (13) 

where f3 is the ratio of radiation pressure acceleration to gravitational acceleration, and 5t 
is the time elapsed since ejection. Particles having a wide range of shapes and compositions are 
approximately represented by /3 = where a^ m is the particle radius expressed in microns 
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(Bohren et al. 1983). At r#= 0.7 AU, the solar gravity is g® = 0.012 m s 2 . Expressing the 


distance 5x in kilometers and the time St in days, we write Equation (13) as 


d/im 


(14) 


In Figure Q, the longest dimension of the debris field on UT 2011 August 28 is of order 100% cor¬ 
responding to Sx ~60,000 km in the plane of the sky. Equation (14) indicates that 1 /im particles 


could have been accelerated over the full length of the debris cloud on timescales as short as St 


1 day. The comet was pictured by Sostero et al. (2011) on Oct. 21 and 23 (St= 64 and 66 days) 
from the ground. A faint fan-like cloud with an extended size l°.5xl0 / .0 was seen in the expected 
position of the comet, corresponding to ~ 10 6 km in the plane of the sky at the geocentric distance 
A ~ 0.25 AU. This extension corresponds, by Equation (14), to the movement of particles up to 
radius a ~ 200 gm in the time since the start of the photometric outburst. Of course, smaller 
particles could have travelled this distance in less time, if the outburst were not impulsive. 

Larger particles, if moving only under the action of radiation pressure, should be confined 
closer to the location of the former nucleus. For example, 10% of the cross-section in the above 
q = 3.3 distribution is carried by particles with a > 200 /im (0.2 mm). Such particles would be 
displaced by radiation pressure over a distance 5x ~ 23,000 km in the NTT image taken on UT 


August 28 ( 5t = 10 days) according to Equation (14). By the same equation, only particles a > 12 
mm would remain within a distance 5x < 350 km 10 days after release, this being the approximate 
size of the seeing disk in the NTT observations of UT 2011 August 28. Such particles carry < 1% 
of the total cross-section, explaining why the distended comet on this date displays no evidence for 
a strong central condensation. 


4.2. Mechanism of the Breakup 

Could the nucleus of Elenin have simply sublimated away? We used Equation ([6]) to estimate 
the thickness of ice that could be sublimated on the pre-perihelion leg of Elenin’s orbit. The 
equation was solved for F s as a function of heliocentric distance, r#, and the equivalent thickness 
of sublimated ice was then calculated from 


A r n 



(15) 


where p n = 1000 kg m~ 3 is the nominal density. We start the integration at time to = 0, corre¬ 
sponding to Elenin at heliocentric distance 5 AU where water ice sublimation is negligible. The 
result of Equations ©> and (|15[) is shown in Figure (^3|, for the two limiting values of parameter 


X- 
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Figure (13) shows that sublimation losses on the journey inwards to perihelion are tiny (roughly 
A r n = 2 to 5 nr, for the two models) compared to the r n ~ 600 nr radius of the nucleus. Therefore, 
sublimation losses cannot be directly responsible for the shrinkage and disappearance of the nucleus 
of Elenin, which requires A r n /r n = 1. This is unlike the case of sun-grazing comets (i.e. those which 
reach perihelion distances of a few solar radii). Many of these bodies have dimensions of order 10 
meters ( |Knight et al. 2010), and are small enough to simply sublimate away (Brown et al. 2011). 
Comet Elenin also disintegrated too far from the Sun for tidal stresses to have played a role (the 
Roche radius for a nucleus of density 1000 kg m -3 is ~ 2.5 Rq ~ 0.01 AU). Gas pressures produced 
by sublimation and ram pressure with the solar wind are likewise both utterly negligible at r# = 
0.7 AU. 

It has been suggested that torques exerted by the loss of material from the nucleus can lead 
to rotational breakup, and that rotational disruption is the likely dominant mode of destruction of 
comets (Jewitt 1992J 1997). Observationally, Bortle (1991) has noted that intrinsically faint (pre¬ 
sumably small nucleus) long-period comets with perihelia <0.5 AU vanish with higher probability 
than bright (presumably large nucleus) long-period comets. Rotational breakup is especially effec¬ 
tive for small cometary nuclei, suggesting that it may be responsible for the demise of Elenin and 
small long-period comets, generally. Accordingly, we examine the hypothesis that the disintegration 
of Elenin was precipitated by a rotational instability. 

Mass lost from a nucleus results in a change in the angular momentum, given by A L = 
kxAMvthi" n , where AM is the mass lost, vth is the outflow speed at which the mass leaves, and r n 
is the nucleus radius. Dimensionless parameter kx is the ratio of the moment arm for the torque 
to the nucleus radius. Simple estimates give kx ~ 0.05 (Jewitt 1997) while measurements from 
two comets give 0.005 < kx < 0.04 (for P/Tempel 1 (Belton et al. 2011)) and kx ~ 0.0004 (for 
103P/Hartley 2 (Drahus et al. 2011)). We take kx ~ 10~ 2±1 , to reflect the dispersion in these values. 
We identify AM with the mass of a spherical shell on the nucleus, AM = 47rr 2 /9 n Ar n , where p n 
is the nucleus density and A r n is the thickness of the shell. Then, AL = \.'Kkxr\p n Vth^ r n- The 
angular momentum of the nucleus is L = koAInr^uj, where fco = 2/5 for a homogeneous sphere of 
mass M n and co = 2 tt/Pq is the angular frequency of rotation at initial period Po- We substitute 
M n = 47rr 3 /9 n /3 and set L = AL to estimate the thickness of the layer which must be lost to modify 
the nucleus angular momentum by a factor of order unity; 


A r r , 


27T fc p r 2 
SkxVthPo 


(16) 


We substitute vth = 500 m s~ 3 corresponding to the sound speed at 200 K, the approximate equilib¬ 
rium temperature of ice freely sublimating at 1 AU. We assume Po = 5 hr, typical of kilometer-sized 
small bodies, kx = 10~ 2±1 and r n = 600 m to find 0.3 < Ar n < 30 m. This is very approximate, 


not least because several of the parameters in Equation (16) are unmeasured. Nevertheless, it is 


permissible to note that for most observationally allowed choices of the parameters, the equation 
gives A r n /r n <C 1, meaning that only a tiny fraction of the nucleus radius needs to be shed in order 
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to change the angular momentum by a large factor. For comparison, the destruction of the nucleus 
by sublimation alone requires A r n /r n = 1, does not lead to a sudden break-up, and takes much 
longer than the time taken by Elenin to travel from its discovery distance to perihelion. 


The range 0.3 < A r n < 30 m is shown in yellow in Figure (13). The Figure shows that both 
the low and high temperature model curves, shown as blue and red lines, respectively, cross into 
the yellow shaded zone at pre-perihelion distances about 1.5 AU and 2.5 AU, respectively. Both 
distances are reached before the date of peak brightness. While this is far from proof that Elenin 
disintegrated because of its own out gassing torques, we must conclude that this is a plausible 
mechanism worthy of further consideration. 
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5. SUMMARY 

We present a study of long-period comet C/2010 XI (Elenin), using observations taken with 
a variety of telescopes over a wide range of heliocentric distances and phase angles. We find that: 

1. The comet displayed a surge in apparent brightness by ~11 magnitudes centered on UT 2011 
August 12.95 (DOY 224.95), when pre-perihelion at heliocentric distance rn= 0.84 AU. The 
peak of this brightness surge is coincident with the passage of the comet through maximum 
phase angle (a =172°.6 on DOY 225), showing that it is caused by strong forward-scattering 
by cometary dust, not by outburst or nucleus disintegration. 

2. A previously undetected, much smaller brightness surge of ~3 magnitudes occurred on UT 
2011 August 30 (DOY 242), about 17 days after the forward-scattering peak. This peak is 
unrelated to the phase angle and instead represents the true breakup of the nucleus, which 
began at heliocentric distance r# = 0.74 AU on UT 2011 August 18±1 (DOY 230±1), some 
23±1 days before perihelion on UT September 10 (DOY 253). The dust cross-section is 
estimated 4300T900 km 2 at the peak of the brightness surge. 

3. The measured gas production rates are consistent with equilibrium sublimation from an area 
~1 km 2 , providing a crude estimate of the pre-disruption radius of the nucleus r n ~ 0.6 km. 
Deep images after the breakup set an upper limit to radius of surviving fragments near 40 m 
(geometric albedo 0.04 assumed). 

4. The increase in the cross-section implied by the brightening is consistent with breakup of a 0.6 
km radius spherical nucleus into a power-law distribution of fragment sizes. We find a power- 
law index q = 3.3T0.2, such that half of the scattering cross-section lies within sub-micron 
dust particles while most of the mass is in large particles, up to the limiting size established 
by CFHT observations. 

5. The pre-perihelion disintegration of Comet Elenin starting at about 0.7 AU cannot be ex¬ 
plained by simple sublimation losses, by gravitational stresses or by sublimation stresses. 
We propose that the nucleus was accelerated to rotational breakup by torques imposed by 
sublimation-driven mass-loss. 


We thank Bernard Hausler and Pedro Lacerda for attempting post-perihelion observations 
on our behalf, and Rachel Stevenson and Jan Kleyna for help with SWARP. We thank Bin Yang, 
Masateru Ishiguro, and Man-To Hui for reading the manuscript and providing comments. We thank 
the referee, Mathew Knight, for his speedy review and valuable comments. We thank the CFHT 
Director for kindly allocating discretionary time at short notice for this work. The Heliospheric Im¬ 
ager instrument was developed by a collaboration that included the University of Birmingham and 
the Rutherford Appleton Laboratory, both in the UK, the Centre Spatial de Liege (CSL), Belgium, 
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Table 1. Observations in 2011 


Dates [UT] 

Telescope 

Field of View 

Pixel Size 

Wavelengths 

Feb 10 - Jun 03 

KAIT 

6'.8 

0.8" 

R-filter 

May 16 - Jun 20 

STEREO HI-1A 

20° 

70" 

630-730 nm 

Aug 01 - 05* 

STEREO HI-2 B 

70° 

240" 

400-1000 nm 

Aug 06 - 12t 

STEREO HI-1 B 

20° 

70" 

630-730 nm 

Aug. 14 - Sep. 04 

STEREO HI-1 B 

20° 

70" 

630-730 nm 

July 23, August 28 

NTT 

4.1' 

0.24" 

R642 nm 

Sep 01 - 22 

STEREO HI-1 A 

20° 

70" 

630-730 nm 

Oct 22 

CFHT 

1° 

0.187" 

R-filter 

Nov 11 

CFHT 

1° 

0.187" 

R-filter 


t During these time periods, STEREO B was specially maneuvered to be rolled toward 
the comet. 
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Table 2. Gas Production Rates and Active Areas 


UT Date 

Q a 

dM/dt h 

r H [AU] c 

F s (min) A 

F s (max) e 


A s {rnax)& 

Reference 

2011 Mar 28 

4.5xl0 26 

12 

2.92 

6.4x 10 -6 

3.0xl0“ 5 

0.4 

1.9 

Korsun et al. (2012) 

2011 Jun 28 

6.3±2.5xl0 27 

190±80 

1.60 

6.6x 10 -5 

1.6x 10 -4 

1.2±0.5 

2.8±1.1 

Schleicher (2015) 

2011 Jul 30 

4.8xl0 27 

144 

1.07 

1.8xl0' 4 

3.9x 10 -4 

0.4 

0.8 

Drahus et al. (2011) 

2011 Sep 07 

<3x 10 27 

<90 

0.49 

9.8xl0~ 4 

2.0xl0- 3 

<0.05 

<0.09 

Lovell et al. (2011) 


a Estimated H2O production rates, molecules s —1 
b Equivalent H2O gas mass production rate, kg s _1 
c Heliocentric distance, AU 

d Minimum equilibrium sublimation flux from Equation § , kg m 2 s 1 
e Maximum equilibrium sublimation flux from Equation |5|, kg m —2 s —1 
f Minimum exposed ice area from Equation (t), km 2 
g Maximum exposed ice area from Equation (|T]l, km 2 
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UT Date (2011) 



Day of Year 


Fig. 1. Heliocentric (black) and geocentric distances (dotted brown, dashed green and dotted- 
dashed blue for Earth, STEREO A and STEREO B, respectively) in 2011. Differently colored 
circles represent observations from individual instruments as labeled in Fig. [2j and also explained 
in Table [U 
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UT Date (2011) 



Day of Year 


Fig. 2.— Phase angle of the comet as viewed from different instruments, each indicated by differ¬ 
ently colored circles. The dotted brown, dashed green and dotted-dashed blue curves represent the 
observations from Earth, STEREO A and STEREO B, respectively. 
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Fig. 3.— Sample image from the STEREO B HI camera taken UT 2011 August 15 at a solar 
elongation of 6.8°. The full image shown is 20° square and has been processed to remove the diffuse 
coronal background but not held stars, so that the density of these can be seen. The inset shows a 
region 23' square centered on comet Elenin. Yellow and green arrows show the projected anti-solar 
direction and the negative of the heliocentric velocity vector. The two brightest point sources are 
the planets Mercury and Jupiter. 
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Fig. 4.— Comet images taken by NTT on UT 2011 July 23 (left) and August 28 (right), showing a 
dramatic physical change. The image panels are 96 // x96 // . Yellow and green arrows show, respec¬ 
tively, the position angles of the projected anti-solar direction and the negative of the projected 
heliocentric velocity vector. North and East are to the top and left, as marked. 
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Fig. 5.— Apparent magnitude of comet Elenin as a function of time expressed as Day of Year (bot¬ 
tom) and UT Date (top). Colored symbols represent the measurements from different instruments, 
KAIT: circles; HI-lBs: diamonds; HI-1B: +; HI-1A: x; and NTT: triangles. The uncertainties are 
estimated from the scatter of repeated measurements at each date. The magnitudes presented in 
this figure and all other figures are R-magnitudes. 
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Fig. 6.— Comet magnitudes after the application of the aperture correction (c.f. Equation [T]) . 
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UT Date (2011) 
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Fig. 7.— Comet magnitudes after the application of both the aperture correction and the inverse 
square law correction (c.f. Equation [ 2 ]) . The peak of the light curve is fitted with a Gaussian profile 
having a center on UT 2011 August 12.95T0.5 (dotted curve). The full width at the half maximum 
of the Gaussian profile is 8 days. 
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Fig. 8.— Same as Figure 0 but plotted versus the phase angle. The pre-peak magnitudes (brown, 
black and blue circles) are fitted by a combined Henyey-Greenstein function (Eq. [4]) to represent 
both forward and backward scattering (solid line). 
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Fig. 9.— Same as Figure ([8]) but separating phase angles before and after the peak. The thick 
black curve is the fitted Henyey-Greenstein function. 
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UT Date (2011) 



Fig. 10.— Absolute magnitudes of Elenin, including aperture correction, inverse square law correc¬ 
tion and phase function correction. The measurements are scaled such that the comet brightness 
at rjj = A=1 and a = 0° is 11.7 magnitude, equal to the best-fit value from the Henyey-Greenstein 
fits. 
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Fig. 11.— Enlargement of Figure (10) near the time of the break-up event. A Gaussian fit gives 
peak brightness on UT 2011 August 30.1T0.5 and the peak brightness reached 2.6 magnitude from 
the base. The full width at half maximum of the Gaussian profile is 13.8 days. 
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Heliocentric Distance [AU] 


Fig. 12.— Sublimating area as a function of heliocentric distance. Blue (red) points correspond 
to models in which the sublimation is assumed to occur at the lowest (highest) equilibrium tem¬ 
peratures for a given heliocentric distance. Arrows signify upper limits based on non-detections of 
gaseous emission lines. 








- 37 



- 7 

- 5 


- 3 


2010.5 


2011.0 2011.5 

UT Date 


2012.0 


X 

o' 

o 

ft 

s 

B. 

o 


1 

0.8 

0.6 

0.4 


a 

o 

re 

70 

> 

C 


Fig. 13.— Integrated thickness of ice layer lost by sublimation on Elenin as a function of the date. 
The red and blue curves correspond to the high and low temperature limits (x = 1 and y = 2, 
respectively) in Equation (15). The yellow shaded region marks the range of ice layer thicknesses 
which must be lost in order to change the nucleus angular momentum by a factor of order unity, 
according to Equation (16). The black short-dashed line shows the heliocentric distance as a 
function of date, for reference (right hand axis). The dates of peak brightness (UT 2011 August 
12) and perihelion (UT 2011 September 10) are marked. 









